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Three main strategies are being pursued to search for non-gravitational dark matter 
signals: direct detection, indirect detection and collider searches. Interestingly, experi¬ 
ments have reached sensitivities in these three search strategies which may allow detec¬ 
tion in the near future. In order to take full benefit of the wealth of experimental data, 
and in order to confirm a possible dark matter signal, it is necessary to specify the nature 
of the dark matter particle and of the mediator to the Standard Model. In this paper, we 
focus on a simplified model where the dark matter particle is a Majorana fermion that 
couples to a light Standard Model fermion via a Yukawa coupling with a scalar mediator. 
We review the observational signatures of this model and we discuss the complementarity 
among the various search strategies, with emphasis in the well motivated scenario where 
the dark matter particles are produced in the early Universe via thermal freeze-out. 
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1. Introduction 

Since Zwicky conjectured the existence of dark mattei!^^ (DM) to explain his own 
observations on the dynamics of galaxies in the Coma Cluster,!^ a lot of effort has 
been invested in determining the properties of this new form of matter. On the one 
hand, current cosmological and astrophysical observations are consistent with the 
dark matter being constituted by an electrically neutral and colorless particle, stable 
on cosmological time scales, non-relativistic at the time of structure formation, and 
with an abundance today which amounts to ^ 80% of the total matter density of our 
Universe. On the other hand, the particle physics properties of the dark matter, such 
as the mass, the spin, the lifetime, or the nature and strength of its couplings to the 
luminous matter, are either unknown or poorly constrained by observations or by 
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theoretical considerations. In fact, viable dark matter models have been constructed 
with masses ranging between ~ 1 /reV and 10^® GeV, and interaction cross sections 
ranging between 10“^® pb and ~ 1 pb (for a review, see Refs. [3||^. 

One of the most pressing open questions in Fundamental Physics consists pre¬ 
cisely in pinning down the Particle Physics properties of the dark matter particle. 
To this end, three different approaches are being pursued: direct detection, which 
aims to the observation of the nuclear recoil induced by the scattering with a dark 
matter particle, indirect detection, which aims to the observation of the photons, 
antimatter particles or neutrinos hypothetically produced in dark matter annihi¬ 
lations or decays, and collider searches, which aim to the production at a collider 
of dark matter particles and their subsequent detection through the observation of 
large amounts of missing energy in the final state. Given the intrinsic complexity of 
searching for a particle with unknown properties in an experiment with large and, in 
some instances, poorly understood backgrounds, it is of paramount importance to 
theoretically identify the cleanest and most powerful search strategy, as well as the 
synergies among the various strategies. This program, unfortunately, relies heavily 
on the details of the model. 

In this review we focus in scenarios where the dark matter particle is a Majorana 
fermion with mass in the range ~ 10 — 10^ GeV that couples to a light Standard 
Model (SM) fermion via a Yukawa coupling with a scalar mediator. This scenario 
leads under reasonable assumptions to the correct dark matter abundance and of¬ 
fers a rich phenomenology both in indirect detection, direct detection and collider 
experiments, which moreover leads to characteristic signals in the concrete scenario 
where the dark matter particle is degenerate in mass with the scalar mediator. In 
Sectionwe present the model and we identify its free parameters. In Sectionwe 
analyze the thermal production of the dark matter particle, emphasizing the role of 
coannihilations in the mass degenerate scenario. In Section]^ we review the signals 
of the model in gamma-ray, antimatter and neutrino detectors, as well as the con¬ 
straints on the model parameters that stem from current searches, while in Sections 
[^andj^the signals in direct detection and collider experiments, respectively. Finally, 
in Section we discuss the complementarity among the various search strategies 
and in Section we present our conclusions. [Appendix A contains a collection of 
relevant cross sections and [Appendix Bj a discussion of effective operators for direct 
detection. 


2. Description of the model 

We consider an extension of the Standard Model by one colorless and electrically 
neutral Majorana fermion, which we assume to constitute the dominant com¬ 
ponent of dark matter in the Universe, and one complex scalar particle, 77 , which 
mediates a Yukawa interaction between the dark matter particle and a Standard 
Model fermion /. The general form of the Lagrangian reads: 

£ n \ r \ r \ /’fermion i /’scalar /i\ 

- >CsM + + ^int • UJ 
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Here, £sm is the Standard Model Lagrangian which includes a potential for the 
Higgs doublet V = On the other hand and are 

the parts of the Lagrangian involving just the Majorana fermion x the scalar 
particle 77, respectively, and which are given by: 




- -X2{v^vf 


( 2 ) 


where denotes the covariant derivative. Lastly, and denote the 

fermionic and scalar interaction terms of the new particles to the Standard Model 
fermions and to the Higgs doublet. The interaction terms depend crucially on the 
quantum numbers of the dark matter particle and the scalar 77. We first impose a Z 2 
discrete symmetry under which x and 77 are odd while the Standard Model particles 
are even, in order to guarantee the stability of the dark matter particle. We also 
impose for simplicity that the dark matter particle only couples to one generation 
of fermions, which can be ensured by postulating that the extra scalar particle 77 
carries a family global quantum number while the dark matter particle does notj^ 
Under these assumptions, the possible interaction terms are restricted by the gauge 
quantum numbers of the dark matter particle and the Standard Model fermion it 
couples to. 

We will focus in this review in scenarios where the dark matter particle is a 
SU{2)l singlet, and therefore does not couple at tree level with the weak gauge 
bosons. In these scenarios, the dark matter hypercharge must be zero in order to 
render an electrically neutral particle. Then, the Yukawa coupling of the dark matter 
to a right-handed quark (lepton) fa with hypercharge Y, requires the scalar 77 to 
be a color triplet (singlet), singlet under SU{2)l and with hypercharge —Y. With 
these quantum numbers the only interaction terms in the Lagrangian are: 


4rr“ = -J/x/i?^ + h.c. 


c 


scalar 


= -A3($'f$)(77'l'77) . 


( 3 ) 


On the other hand, if the dark matter only couples to a left-handed quark (lepton) 
with hypercharge Y, the scalar 77 must be a color triplet (singlet), doublet under 
SU{2)l and with hypercharge —Y. The corresponding interaction terms are: 


C- 


fermion 


c 


scalar 


= -yx{fLi<J 2 V) + h.c. , 


( 4 ) 


In the Minimal Supersymmetric Standard Model (MSSM) these possibilities are 
realized if x is the bino and 77 is a sfermion. The Yukawa coupling y = t/mssm is 
then fixed in terms of the U{\)y gauge coupling g' to t/mssm = 4(7V(3'\/2) ^ 0.33 


^Lifting this requirement leads in general to too large flavor violating effects. A detailed discussion 
of the conditions that must be fulfilled in order to satisfy the constraints from flavor physics can 
be found in the Appendix of Ref. for coupling to quarks, and in Ref. I^for coupling to leptons. 
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for up-type right-handed quarks, j/mssm = ^ 0.16 for down-type right- 

handed quarks and j/mssm = g '^ 0.08 for the left-handed quarks. Besides, 
2 /mssm = V^g' ~ 0.48 for right-handed leptons and yussM = g'lV^ ~ 0.24 for 
left-handed leptons. For a general mixed neutralino, the couplings to the left- or 
right-handed fermions are, respectively, given bjlH 

2/mssm,l = — g'Nii =F -^Ni2 , (5) 

2/mssm,r = Aiqfg'Nil ■ (6) 


In the following we explore the parameter space of the simplified model treating 
the coupling j/ as a free parameter. Nevertheless, we refer to the supersymmetric 
coupling strength as a benchmark value. Furthermore, for simplicity, we assume 
that I Ail <C 1, such that the relevant parameters of the model are the dark matter 
mass the mediator mass rriri and the Yukawa coupling y. In addition, we mostly 
focus on a coupling to the right-handed fermions for definiteness, although many 
results can be generalized in a straightforward manner. 

Apart from the collider constraints discussed in detail below, the Yukawa in¬ 
teraction involving leptons gives rise to a contribution to the anomalous magnetic 
moment 


Nat = - 


y 


„2m® 
-m|-A. 








(7) 


IBtt^ 6(m2 — 

For both electrons and muons the contribution is subdominant compared to the 
Standard Model expectation for Yukawa couplings compatible with indirect detec¬ 
tion constraints and > 10 GeVim 


3. Relic abundance 


One of the most attractive frameworks to generate the observed dark matter abun¬ 
dance in our Universe is the freeze-out mechanism. It assumes that the dark matter 
has interactions with the Standard Model strong enough to keep at high temper¬ 
atures the dark matter particles in thermal equilibrium with the Standard Model 
particles. Besides, the interactions must be weak enough to allow the dark mat¬ 
ter particles to drop out of equilibrium at sufficiently early times. After this time, 
dubbed the freeze-out time, the number density per comoving volume of dark mat¬ 
ter particles remains practically constant until today, thus constituting the dark 
matter population we observe in the Universe. Under these assumptions, the dark 
matter abundance = 0.1198 ± 0.0015,^ as measured by the Planck satellite, 

can be naturally explained. 

Defining the variable cct.o. = temperature Tf.o. at which the freeze- 

out takes place is approximately given by the solution of the following equation®! 


.O. 


log 


0.038c(c -|- 2)gintMpim^{av) 


g*{xf.oy^'^Xi' 


1/2 


( 8 ) 
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where c is a numerical factor, c « 0.5, Mp\ = 1.22 x 10^® GeV is the Planck mass, 
g* ixi,o .) is the number of relativistic degrees of freedom at the freeze-out tempera¬ 
ture and gint is the number of dark matter internal degrees of freedom, in this case 
gint = 2. Typically Xf.o, =20-30. 

The relic density of dark matter particles can be calculated solving the Boltz¬ 
mann equation for the number density of dark matter particles as a function of the 
temperature. The result approximately readsl^^ 


2_ 1.07 X lO^GeV^ 


(9) 


where T(a;f,o.) = {av)x~‘^ dec . Here {av) is the thermally-averaged annihilation 

cross section, which can be readily calculated using the formalism presented in 
Ref. When expanding the cross section for low velocities as av = a + the 
thermal average is given by (av) = a + 3b/x{,o.- 

As apparent from the previous equations, the thermally-averaged annihilation 
cross section is crucial for determining the dark matter relic abundance. For a 
SU{2)L singlet dark matter particle, the dominant annihilation channel is XX //, 
via the t-channel exchange of the mediator 77 . The velocity independent part is given 

by 


a = 


Nc mj - mj 

32 TT — mj + rn'fj 


2 > 


( 10 ) 


and vanishes in the limit to/ —)■ 0. Besides the velocity dependent part approxi¬ 
mately reads, 


Ncm'^y'^ (to^ -I- 
487r (to^ -I- TO^) 4 


( 11 ) 


in the limit to/ —>■ 0 (see Appendix A| for the general expression). 

For light fermions, namely to/ ^ to^, and to^/to;^ > 1.5 the relic density is 
approximately given by 




0.12 

70.37^ 

( 

V 

1 - 

X 

-f 

1—1 

Nc 

1 y / 

aoo GeV/ 

. i 

(1 - 1 - 


( 12 ) 


for a freeze-out temperature Tf.o, = to ^/20 and g*{T{,o.) = 80.111^ 

This formalism must be modified when the scalar mediator is very degenerate 
in mass with the dark matter particle. If this is the case, the scalar y is still present 
in the thermal bath and depletes the number of dark matter particles. The effective 
cross-section is given by 


(XeflV 


eq eq 

> , 7 ^^—^777x0-,, 

ij 




eq \2 ’ 


(13) 
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where the sum extends over all coannihilating particles, including x- Here, aijV is 
the cross section for each of the reactions participating in the freeze-out, shown 
in Fig. 0 and = gi(TOiT/(27r))^/^e is the number density per comoving 

volume of particles i with mass rrii in equilibrium at the temperature T, with = 2, 
9y = 9fj = 3(1) being the number of internal degrees of freedom of the dark matter 
and of the colored (uncolored) scalar mediator. 

The requirement of reproducing a dark matter density with the value measured 
by the Planck satellite, = 0.1198±0.00ld^ fixes one of the parameters of the 

model in terms of the remaining two, e.g. the Yukawa coupling y = 

Under this reasonable condition, the parameter space of the model is spanned by 
only two parameters, that can be chosen to be rriy^, the dark matter mass, and 
{rriri — which measures the mass difference between the scalar mediator 

and the dark matter particle. Fig. shows contour lines with the value of the 
Yukawa coupling that leads to the observed dark matter abundance via thermal 
freeze-out, assuming that the dark matter particle couples to leptons (left plot) 
or to quarks (right plot). As apparent from the plot, the region with large 
and/or large leads, for perturbative couplings y < -s/Ur, to a relic abundance 
which exceeds the observed value. On the other hand, the region with small 
and/or large degeneracy leads to a too small relic abundance, due to the crucial 
role of coannihilations during freeze-out. This behavior can be understood from the 
expression of the effective thermal cross section in the coannihilation regime, Eq.[T3| 
In this simplified scenario, it is approximately given by: 

aeffV ~ R + (Trir^V , (14) 


which includes a Boltzmann suppression factor R = n®'^/n®'^ oc e "’i’ for the 
processes involving the coannihilating particle 77. The cross sections for %%, xi! 
and rjr] are proportional to j/^, y'^g'^ and g'^, respectively, g being a gauge coupling 
constant (see Fig.[^. Besides, for fixed ratio the cross sections must be 

inversely proportional to m^. Then, casting the cross sections as 

one findJi^ 


XI (X ^ ^ _ flSl 

^ (TTeffu) yHC^^)+y^g^C^y)+gHCr,r,) ' ^ ^ 

As apparent from this equation, there is always a value of the dark matter mass 
below which the dark matter abundance lags below the observed value, even for 
y = 0, provided g^Cr^rj 7 ^ 0 , which is the case in the coannihilation regime where the 
Boltzmann suppression is relatively mild. In very degenerate scenarios, m^/m^ « 
1.1, the lower limit on the dark matter mass is > 200 GeV (50 GeV) for dark 
matter coupling to quarks (leptons). 
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Fig. 1. Feynman diagrams relevant for the computation of the relic density, including coannihila¬ 
tions with a colored mediator. For an uncolored mediator, similar diagrams involving SU{2) x U{1) 
gauge bosons are relevant. 


DM coupling to leptons DM coupling to quarks 



[GeV] m, [GeV] 


Fig. 2. Parameter space compatible with dark matter production via thermal freeze-out for 
coupling to leptons (left plot) and to quarks (right plot). The contour lines show the value of the 
Yukawa coupling y required to reproduce the observed dark matter density = 0.1198 di 

O.OOlP via thermal freeze-out. The upper right shaded region leads to > Hom unless the 
coupling becomes non-perturbative, y > \/47r. On the other hand, the lower left shaded region 
leads to < floM due to strong coannihilations. The plot also shows the lines corresponding 
to the MSSM values of the Yukawa coupling, namely y = 0.33 (0.5) for a right-handed squark 
(slepton) mediator. 


4. Indirect detection 

Immediately after the dark matter freeze-out, the annihilation rate of dark matter 
particles became negligibly small and their number density per comoving volume re¬ 
mained practically constant. In this epoch, the dark matter distribution was roughly 
homogeneous and isotropic throughout the whole Universe, except for the small den¬ 
sity fluctuations generated by the inflaton field. The overdensity regions accreted 
more and more dark matter particles, thereby generating regions in our Universe 
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with a density of dark matter particles which became orders of magnitude larger 
than the average density. The baryonic matter followed the gravitational potential 
wells generated by the dark matter, leading to structures which are identified with 
the galaxies and clusters of galaxies we observe today. In these regions of enhanced 
dark matter density, annihilations may still be efhcient, producing Standard Model 
particles which contribute to the total fluxes of gamma-rays, cosmic-rays and (anti- 
)neutrinos. The search for this exotic component over the expected astrophysical 
backgrounds then provides information about the dark matter properties. 

The annihilation of dark matter particles with mass produces Standard 
Model particles at the position r, at a rate per unit kinetic energy T and unit 
volume given by 


Q(T,rl 


1 

2 


'^icrv)i 


dN^ 

It 


(16) 


where {av)i is the annihilation cross section times the relative dark matter velocity 
in the annihilation channel i, producing an energy spectrum of Standard Model par¬ 
ticles dN'‘/dT, and is the dark matter density at the position r. Unfortunately, 
in this expression, none of the parameters involved is positively known. 

To assess the prospects to observe a dark matter signal in the sky, it is common 
to fix the annihilation final state, which then fixes dN/dT. Besides, the distribu¬ 
tion of dark matter particles in the Milky Way is not precisely known, although it 
can be inferred from numerical iV-body simulations. Some popular choices for the 
dark matter density profile, which illustrate the range of uncertainty in the predic¬ 
tions of the indirect dark matter signatures, are the Navarro-Frenk-White (NFW) 
profile Enin! 


Pdm(»’) = 


Po 


(r/rs)[l-b (r/rs)]2 ’ 

with scale radius rg = 24 kpc,li^ the Einasto profileE^HIll 

^2 


(17) 


/3dm(?’) = Poexp 


a 


(18) 


with a = 0.17 and Vg = 20 kpc,^and the much shallower isothermal profile?^ 


PDuir) = 


Po 


1 -f r^/r^ ’ 


(19) 


with Tg = 4.4 kpc. In all the cases, the overall normalization factor po is chosen to 
reproduce the local dark matter density p© = 0.39 GeV /cm^^oHlll at the distance 
r© = 8.5 kpc of the Sun to the Galactic center. 

Under these assumptions for the annihilation final state and the dark matter 
density distribution, the source term Eq. (16) depends on two free parameters, the 
dark matter mass and the annihilation cross section in that channel {(tv), which 
span the parameter space of indirect dark matter searches. Various experiments 
are currently constraining this parameter space using gamma-rays, antimatter or 
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neutrinos. Let us first discuss the possible annihilation channels of the model and 
then the methods for indirect searches in each of these channels, as well as the 
present experimental limits. 


4.1. Annihilation channels 


As discussed in the previous section, the cross section for the lowest order anni¬ 
hilation channel XX f f is suppressed by the mass of the final fermion and by 
the relative velocity of the dark matter particles in the Galactic center, v « 10“^. 
Therefore, higher order processes that lift the helicity suppression can give a sizable 
or even dominant contribution to the total annihilation cross section for relative 
velocities relevant to indirect detection. This is the case for the two-to-three annihi¬ 
lation into a fermion-antifermion pair with the associated emission of a gauge boson 
XX //vESEH and the one loop annihilation into two gauge bosons, xx —>■ VV, 
which have the largest branching fractions over most of the parameter space, de¬ 
spite the suppression of the cross section by the additional gauge coupling and the 
three-body phase space for the former process, and by the loop factor for the lat¬ 
ter. The two-to-three annihilation xx ffh^ with h the Higgs boson, also lifts 
the helicity suppression, however they are subdominant unless the mediator couples 
strongly to the Higgs, Xi ~ 0{1) (notice that the s-wave amplitude for the one-loop 
annihilation xx hh, contrary to xX W, identically vanishes due to CP con¬ 
servation). Explicit expressions for the annihilation cross sections for the various 
channels can be found in Appendix A| 

The relative importance of the different final states mostly depends on the masses 
of the dark matter particle and the scalar mediator. This is illustrated in Fig. 
which shows the dependence of the branching ratios with the dark matter mass 
(upper plots) and with the mass degeneracy parameter mrj/m^ (lower plots), 
for dark matter coupling to cr (left plots) and to ur (right plots), assuming that 
all quartic couplings vanish. As apparent from the upper plots, sufficiently away 
from the production threshold the branching ratios for the two-to-three processes 
and the one loop processes have constant branching ratios; only the two-to-two tree 
level annihilation into // shows a strong dependence on the mass. Moreover, it 
follows from the lower plots that when the mass difference between the scalar medi¬ 
ator and the dark matter is very large, the two-to-three processes become negligible 
and the two-to-two dominate the annihilation. Among these, the process with the 
largest branching fraction is the annihilation into 77 , for coupling to leptons, and 
into gg^ for coupling to quarks. In contrast, as the scalar mediator becomes more 
and more degenerate in mass with the dark matter particle, the branching ratios 
for the two-to-three processes become larger and larger, and even become the dom¬ 
inant annihilation channels in the very degenerate limit. The largest cross section 
in this regime corresponds to the annihilation into //q, for coupling to leptons, 
and into ffg, for coupling to quarks. The ratios of strong and electroweak internal 
bremsstrahlung cross sections relative to //q, as well as gg relative to 77 , are shown 
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Branching fractions for coupling to eji and = 2 Branching fractions for cou])ling to ur, and rurj/m,^ = 2 




Branching fractions for coupling to and = lOOOGeV Branching fractions for coupling to ur and = lOOOGeV 




Fig. 3. Branching ratios of the various tree-level and one-loop annihilation channels for 
dark matter coupling to the right-handed electron (left plots) or up-quark (right plots). 
The top plots show the branching ratios as a function of the dark matter mass for fixed 
mrjlm^ = 2, while the bottom plots, as a function of rrir^lm^ for fixed = 1000 GeV. 


in Tablethese values are valid for m/ —0 and :$> mzl‘^- 


Table 1. Ratios of annihilation cross sections for a scalar mediator 
that couples to right-handed up-type quarks, down-type quarks and 
leptons, respectively, in the limit my ^ 0. The ratios for XX ff^ 
are valid for rriy^ ^ mz 1"^- In the numerical values the strong coupling 
is evaluated at a scale /i = 300 GeV. 





dR 




38.4 

12as ^ 

154 


OCem 




^^ffz 

tsin^{6w) 

~ 0.30 

tsin^{6w) 

~ 0.30 tan2(e„) ~ 0.30 



; 184 

ISal 

2944 
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4.2. Gamma-rays 


Gamma-rays from dark matter annihilations propagate through our Galaxy in 
straight lines and practically without losing energy. This allows to devise several 
search strategies which exploit the spatial distribution expected from the signal 
events, and which is not possible in antimatter searches. In particular, the most 
prominent targets for indirect searches for dark matter annihilations with gamma- 
rays are the Milky Way center and dwarf spheroidal galaxies. Furthermore, the fairly 
good energy resolution of gamma-ray telescopes allows to search for the character¬ 
istic gamma-ray spectrum produced in a given annihilation channel, an approach 
which is at present not possible with neutrino searches. For these reasons, gamma- 
rays from self-annihilations are a promising and unique channel for dark matter 
detection. 

As discussed in the previous section, in the limit of zero relative velocity, the 
annihilation into a fermion-antifermion pair is helicity- and velocity- suppressed, 
hence the dominant annihilation channels are the two-to-three process XX f 
and the one loop process yy —)■ VV, with V,V' gauge bosons (see Fig. |^. Of 
particular importance for indirect searches are the final states with V = x, since they 
produce a sharp feature in the gamma ray spectrumB^B^IMtl2l More specifically, 
the two-body annihilation yy —)■ 77 produces two monoenergetic photons at = 
m^, hence the energy spectrum is markedly different to the smooth astrophysical 
background. A related process is the annihilation yy —)■ 7 Z, which produces one 
monoenergetic photon at Ej = m^(l —/4 to^). The gamma-ray line generated in 
the 7 Z final state is fainter than the line from 77 , due to the smaller branching ratio 
of this process and the smaller multiplicity of photons in the final state, and can be 
neglected in most instances within the simplified models considered in this review. 
Finally, when the mass difference between the dark matter particle and the scalar 
mediator is small, the two-to-three annihilation yy —)■ //y produces a sharp feature 
in the spectrumpHHH due to the enhancement of the scalar t-channel propagator in 
the non-relativistic limit of the incoming dark matter particles when the momentum 
of the final fermion is close to zero, and which corresponds to an energy of the photon 
close to the kinematical end-point of the spectrum of annihilations. Under the above- 
mentioned conditions, this process generates a signal that resembles a distorted 
gamma-ray line and cannot be distinguished from the spectrum of monochromatic 
photons with the energy resolution of current gamma-ray telescopes. Many works 
have recently investigated the signatures from internal bremsstrahlung, see e.g., 
Refs. [8l [T^[25l|3lH38| In either case, the observation of a sharp feature in gamma- 
rays would constitute a strong hint for dark matter annihilations as no astrophysical 


process can mimic such a signal (a possible exception has been proposed in Ref. 391. 

The total gamma ray spectrum consists of three contributions: the monochro¬ 
matic gamma-ray lines from two-body annihilation with one or two photons in the 
final state, the sharp spectral feature from electromagnetic internal bremsstrahlung, 
and the soft spectrum of secondary gamma rays produced in other annihilation 
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Fig. 4. Feynman diagrams for the annihilation channels which lead to pronounced spectral 
features close to the end-point E^/ = m^. The first three contribute to the process of internal 
bremsstrahlung, while the last two, to the generation of gamma-ray lines. 


channels. Some exemplary spectra are shown in Fig.[^ for the cases nir^/m^ = 1.01 
(upper left plot), 1.1 (upper right plot), 2 (lower left plot) and 5 (lower right plot), 
assuming coupling to the right-handed up-quark and an energy resolution of 10 %, 
which is typical for current gamma-ray telescopes. In the plot, we have neglected 
the contribution from annihilations into ■jZ which, as mentioned above, gives a sub¬ 
dominant contribution to the total gamma-ray flux. It is noticeable from the figure 
the change in the relative weight of the internal bremsstrahlung processes and the 
one-loop processes in the photon spectrum, and which follows from the dependence 
of the branching ratios on cf. Fig. Namely, in the degenerate limit the 

internal bremsstrahlung dominates the total spectrum of annihilations at the high¬ 
est energies, while in the hierarchical case, it is the one loop annihilation into 77 . 
Besides, the relative importance of the annihilations into uug and unj on one hand, 
and gg and 77 on the other hand, is fixed (see Fig. and Table [^. For coupling to 
down-type quarks, the feature from internal bremsstrahlung is suppressed compared 
to the continuum by a factor of 4, and the monochromatic contribution by a factor 
of 16, due to the smaller electric charge of the down-quark. For coupling to leptons, 
on the other hand, the contribution from continuum gamma rays is much smaller 
and the sharp spectral features are even more salient. 

The gamma-ray flux produced in dark matter annihilations is dominated close 
to the kinematical endpoint by the channels XX ffl XX 77 (^-s mentioned 
above, we neglect in this discussion the contribution from 7 Z, which is suppressed 


with respect to 77 ). More specifically, the source term Eq. 16 reads in this case: 




d{c 


Jgqi 


dE^ 


-\- 2{av)^^d{Ej — m^) 


( 20 ) 


Hence, the differential gamma-ray flux as seen from the Earth under the angle ^ 
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Fig. 5. Gamma ray spectrum from dark matter annihilation for a colored mediator coupling to 
right-handed up-quarks, for = 1.01 (left upper plot), rriy^lmy^ = 1.1 (right upper plot), 

rrirjlm^ = 2 (left lower plot), rrinlmy, = 5 (right lower plot). 


with respect to the Galactic Center is given by, 


dEydfl 


1 

Itt 


ds Q{Ej, r) 


( 21 ) 


where r = \/{tq — s cos + (* sin and ro = 8.5 kpc. 

Searches for line-like features have been conducted with data collected by the 
Fermi-LATHEni in the energy range 5-300 GeV and by the H.E.S.S. collaboratioii^ 
in the energy range 500 GeV-25 TeV. The 95% C.L. upper limits on the anni¬ 
hilation cross section + 2av^~j obtained from the combined contribution of 

internal bremsstrahlung and lines were calculated in Ref. and are shown as a blue 
line in Fig.[^ upper plots, for dark matter particles coupling to the right-handed 
muons (left plot) or to the right-handed up quarks (right plot), and for a mass ratio 
with the scalar mediator m^/m^ = 1.01, 1.1 or 2. The limits range from about 
5 X 10“^®cm^/s to 10“^®cm^/s for dark matter masses from 40 GeV up to 10 TeV, 
and depend relatively weakly on the mass ratio The plot also shows the 

corresponding value of the cross section into sharp gamma-ray spectral features ex- 
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pected for a thernially produced dark matter particle, which lags well below the 
current limits, although for some choices of parameters it is only necessary a boost 
in the flux from annihilations by a factor of 5-10 to produce an observable sig¬ 
nal. The next-generation of gamma ray telescopes GAMMA-40CP^ and especially 
CtA,S^ will improve the limits by a factor of a few, depending on the energy, as 
shown in Fig.[^ lower plots. 

We also show in the figures for comparison the present limit and future reach in 
the annihilation cross section into sharp spectral features that can be derived from 
the limits on the model parameters from current and projected direct detection 
experiments, which depend strongly on the mass difference between the mediator 
and the dark matter particle. The derivation of these limits will be discussed in 
more detail in Sec.[^ It is notable from the plot the strong impact of direct detec¬ 
tion experiments on the searches for sharp spectral features, especially in the mass 
degenerate scenario when the dark matter particle couples to the right-handed up 
quark. 

Finally, complementary limits on the model can be derived from the search for 
the continuum emission of gamma-rays in the annihilation. Among the existing 
targets, dwarf spheroidal galaxies are particularly suitable for the search, due to 
their very large mass-to-light ratio and their moderate astrophysical activity, which 
translates into a negligible gamma-ray emission from astrophysical sources. Current 
data show no significant excess of gamma-rays in the direction of dwarf galaxies, 
thus allowing to set fairly stringent upper limits on the flux.®llZl The flux limits 
can in turn be translated into limits on the annihilation cross section in the channel 
XX f fV and which can produce a continuum flux of gamma-rays. These limits 
can be particularly relevant when the dark matter particle couples to a quark, due 
to the large branching fraction of the annihilation XX QQ9- In case, the 
continuum flux of gamma-rays is generated in the fragmentation and decay of the 
gluon and the quarks. The impact of the dwarf limits on the model is illustrated in 
Fig.|^ upper left plot, for the particular case q = ur (the limits for q = da are a 
factor of four stronger, due to the different quark charges). Here, the dwarf limits on 
the model parameters have been translated into limits on the cross section for the 
annihilations producing sharp spectral features. As apparent from the plot, for a 
colored mediator dwarf galaxy observations provide limits on the model parameters 
which can be competitive to the limits obtained from the search for sharp spectral 
features. For an uncolored mediator, on the other hand, the continuum part of the 
annihilation spectrum is much fainter due to the absence of gluon bremsstrahlung, 
therefore, in this case the strongest limits stem from the search for sharp spectral 
features. 


4.3. Antimatter 

Dark matter annihilations in general produce antimatter particles that could be 
detected at the Earth as an excess over the expected astrophysical backgrounds. 
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DM coupling to RH muon (limits) 



DM coupling to RH muon (prospects) 



DM coupling to RH up-quark (limits) 



DM coupling to RH up-quark (prospects) 



Fig. 6. Current constraints from Fermi LAT and H.E.S.S. (upper row) and future 
prospects for CTA and GAMMA-400 (lower row) on the spectral gamma ray feature 
resulting from the superposition of internal bremsstrahlung and loop-induced annihilation 
into monochromatic photons, for mr^/m^ = 1-01 (solid), 1.1 (dashed) and 2 (dotted). The 
left column corresponds to a charged mediator which couples to leptons (here and the 
right column to a colored mediator which couples to quarks (here ur). For comparison we 
also show the cross section expected for a thermal relic for each value of mrj/m^ (black), 
the cross-section for right-handed squark or slepton mediators (gray), as well as current 
constraints from XENONIOO/LUX and prospects for XENONIT (red, see Sec.j^. 


constituted by a secondary component from spallation of high energy cosmic rays 
on the interstellar medium, and possibly a primary component from sources, such as 
the positrons produced by the interactions of high-energy photons in the strong mag¬ 
netic fields of pulsar^^^M^ or the positrons and antiprotons produced by hadronic 
interactions inside the same sources that accelerate galactic cosmic raysEUlUl 
Antimatter particles have electric charge and propagate through the Milky Way 
in a complicated way before reaching the Earth. The propagation of charged parti¬ 
cles in the Galaxy is commonly modeled by a stationary two-zone diffusion model 
with cylindrical boundary conditions. In this scheme, the number density of an¬ 
tiparticles at the position r per unit kinetic energy, f(T,r,t), satisfies the following 
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Table 2. Astrophysical parameters compatible with the 
B/C ratio that yield the minimal (MIN), median (MED) and 
maximal (MAX) antiproton fluxes from dark matter annihi¬ 
lations; taken from Ref. |54| 


Model 

<5 

Kq (kpc^/Myr) 

L (kpc) 

Vc (km/s) 

MIN 

0.85 

0.0016 

1 

13.5 

MED 

0.70 

0.0112 

4 

12 

MAX 

0.46 

0.0765 

15 

5 


transport equation!^ 

0 = = V • [K{T, r)Vf] + ^ [biT, r)f] - V • [V,{r)f] - 2M(z)ran„/ + Q{T, r) . 

( 22 ) 

The boundary conditions require the solution f{T,f^t) to vanish at the boundary 
of the diffusion zone, which is approximated by a cylinder with half-height L = 
1 — 15 kpc and radius R = 20 kpc. 

The first term on the right-hand side of the transport equation is the diffusion 
term, which accounts for the propagation through the tangled Galactic magnetic 
field. The diffusion coefficient K{T,r) is assumed to be constant throughout the 
diffusion zone and is parametrized by K{T) = Kq (3 TZ^, where P = v/c and 
TZ is the rigidity of the particle, which is defined as the momentum in GeV per 
unit charge, TZ = p{GeV)/Z. The second term accounts for energy losses due to 
inverse Gompton scattering (IGS) on the interstellar radiation field (ISRF) as well 
as synchrotron radiation and ionization. The third term is the convection term, 
which accounts for the drift of charged particles away from the disk induced by 
the Milky Way’s Galactic wind. It has axial direction and is also assumed to be 
constant inside the diffusion region: I4(r) = Vc sign(z) k. The fourth term accounts 
for antimatter annihilations with rate Tann when it interacts with ordinary matter 
in the Galactic disk, which is assumed to be an infinitely thin disk with half-width 
h = 100 pc. Lastly, Q{T,r) is the source term of antimatter particles, given in 
Eq. |16| In this equation, reacceleration effects and non-annihilating interactions of 
antimatter in the Galactic disk are neglected. 

The transport equation, using the parametrizations of the different terms given 
above, has a number of free parameters which are inferred from measurements of flux 
ratios of primary and secondary cosmic-ray species, mainly the Boron-to-Carbon 
(B/G) ratio. Unfortunately, the model parameters cannot be determined uniquely 
from current observations, due to degeneracies among the parameters entering in the 
calculation of the cosmic-ray fluxes, for instance between the diffusion coefficient 
and the height of the magnetic diffusion zone. The ranges of the astrophysical 
parameters that are consistent with the B/G ratio and that produce the minimal 
(MIN), median (MED) and maximal (MAX) antiproton fluxes from dark matter 
annihilations were calculated in Ref. [54] and are listed in Table |2l 
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Finally, the flux of primary antiparticles at the Solar System from dark matter 
annihilations is given by: 




(23) 


where v is the velocity of the antimatter particle. 

At energies smaller than ~ 10 GeV the antimatter fluxes at the top of the 
Earth’s atmosphere can differ considerably from the interstellar fluxes due to so¬ 
lar modulation effects. One frequently used parametrization of the effect of solar 
modulation, which can be derived from the full diffusion and convection equations 
describing the solar wind, is the force-field approximationHHUll The fluxes at the 
top of the atmosphere in this approximation are related to the interstellar fluxes by 
the following relation!^ 




(24) 


where m is the mass of the cosmic-ray antimatter particle and Tjs = Txoa + '(’F; 
with Tig and TpoA being the kinetic energies of the antimatter particles at the 
heliospheric boundary and at the top of the Earth’s atmosphere, respectively, and 
(j>F being the Fisk potential, which varies between 500 MV and f.S GV over the 
eleven-year solar cycle. 

For antiprotons, due to their comparatively large mass, energy losses are negli¬ 


gible and therefore the general transport equation, Eq. (221, can be simplified. The 


transport equation for the antiproton density, fp{Tp,r,t) then reads 

0 = ^ = V • {K{Tp, r)Vfp) - V • (y,{r)fp) - 2h5{z)T,^^fp + Q{Tp, r) , (25) 


where the annihilation rate, Fann, is given by 


Tann = («H + 4^/^nHe)o'|p"r)p . 


(26) 


In this expression it has been assumed that the annihilation cross-section between an 
antiproton and a helium nucleus is related to the annihilation cross-section between 
an antiproton and a proton by the simple geometrical factor 4^/^. Furthermore, 
riH ~ 1 cm“^ is the number density of Hydrogen nuclei in the Milky Way disk, nue ~ 
0.07 riH the number density of Helium nuclei and is the proton-antiproton 
annihilation cross-section, which is given e.g. in Refs. [58l[5^ 

The spectrum of cosmic antiprotons has been measured by WiZ- 
ard/CAPRIGE,l^ between 0.62 and 3.19 GeV, AM^ between 0.2 and 4 GeV, 
and more recently by the PAMELA satellit^^ between 60 MeV and 180 GeV and 
BESS-Polar LP^ between 0.17 and 3.5 GeV. In the near future, data from the AMS- 
02 experiment are expected to yield further information on the cosmic antiproton 
fluxEMU] The measured flux as well as the antiproton-to-proton ratio agree well 
with the expectations from secondary production of antiprotons from spallation of 
cosmic ray nuclei, mainly protons and Helium, on the interstellar mediumUBEZl This 
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allows to set stringent upper limits on a possible primary contribution generated 
from dark matter annihilations. 

Annihilations of Majorana dark matter particles coupling to leptons do not pro¬ 
duce antiprotons in the lowest order process XX ^^£“,with £ = e, t. Neverthe¬ 
less, the higher order process xX t~ Z produces antiprotons in the fragmenta¬ 
tion of the Z boson, and therefore the antiproton measurements can constrain this 
model, despite its leptophilic nature. This process has been studied in various works, 
e.g. Refs. Fig. shows the 95% C.L. upper limits on the annihilation cross 

section for the higher order process xX Z calculated in Ref.|A^from the non¬ 


observation of a significant excess of antiprotons in the PAMELA p/p datsP^ over 
the background flux calculated in Ref. In the plot it was assumed the Einasto 
profile and mri/m^ = 1-1, however the limits are rather insensitive to the choice of 
halo profile and to the concrete value of the mass of the scalar mediator as long 
as ~ 0[X). Besides, in order to bracket the astrophysical uncertainties, the 

plot shows limits for the three propagation models discussed above. For the MED 
propagation model the limits are rather weak and lie in the range av < 10“^® cm^/s 
for ra^ = 100 GeV and 10“^^ cm^/s for = 1 TeV, which approximately translate 
into an upper limit on the Yukawa coupling y < 5.4 for = 100 GeV and y < 26 
for = 1 TeV; these limits are very sensitive to the choice of the propagation 
parameters, as apparent from the figure. The plot also shows for comparison the 
limits on that same channel inferred from the non-observation of a significant ex¬ 
cess of gamma-rays in the Fermi-LAT® and H.E.S.S.I^ data, which are independent 
of the propagation model. It follows from the plot that for leptophilic models the 
limits on the cross section derived from gamma-ray observations are more than one 
order of magnitude stronger than those derived from antiprotons. 

On the other hand, the lower left and right plots show, respectively, the corre¬ 
sponding limits for the annihilation channel xX xX ddg. The limits 

on the cross section are only a factor of ^ 3 better than for yy —>• £~^£~Z, however 
the limits on the Yukawa coupling y are significantly better for dark matter parti¬ 
cles coupling to quarks, due to the larger cross section of the three body final state 
involving one gluon. More specifically, for the MED propagation model, it follows 
that y < 0.8 for = 100 GeV and y < 5.4 for = 1 TeV. The plot also shows 
the limits inferred from the Fermi-LAT and H.E.S.S. data, and which are a factor of 
four weaker for couplings to down-quarks compared to the couplings to up-quarks 
due to the different electric charges. As apparent from the plot, for couplings to 
quarks the limits from antiprotons are competitive, and in some instances better, 
than the limits from the search for gamma-ray spectral features. 

The production of antiprotons and antineutrons in dark matter annihilations 
leads to the possibility of also producing antinuclei, such as antideuteron^ or 
antiheliuml5iJl®21 These channels are particularly interesting due to the very small 
background flux of antideuteronJUUHl and antiheliurrEUUll expected at the energies 
relevant for experiments. In fact, the expected flux from spallations of cosmic rays 
on the interstellar medium lies more than three orders of magnitude below the 
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DM coupling to leptons 



DM coupling to u-quarks DM coupling to rf-quarks 




Fig. 7. 95% C.L. upper bounds on the annihilation cross section for the higher order processes 

XX —t l '^l~Z (upper plot) XX uug (lower left plot) and XX ddfl (lower right plot) derived in 
Ref. [^froni the PAMELA data on the antiproton-to-proton fraction|®l adopting the MIN, MED 
and MAX propagation models defined in Table|^ In the figures, it was assumed the Einasto profile 
and a dark matter particle quasi-degenerate in mass with the scalar mediator. For comparison, 
we also show as dashed lines the limits on the cross sections for the corresponding channels from 
gamma-ray data from the Fermi-LAT® and H.E.S.S.EH data. 


current limit by BES^^ and more than one order of magnitude below the expected 
sensitivity of AMS-02 and GAPS,l^therefore, the observation of a few antideuterons 
in experiments will constitute a strong hint for their exotic origin. The antideuteron 
flux at the Earth from dark matter annihilations is, on the other hand, strongly 
constrained by the non-observation of an excess in antiprotons.^^ Given the stringent 
limits existing from PAMELA, and despite the various sources of uncertainty in the 
modeling of the antideuteron production, the observation of an antideuteron flux 
at AMS-02 or GAPS from dark matter annihilations will be challenging.!^ 

Dark matter annihilations also produce positrons. In this case galactic convection 
and annihilations in the disk can be neglected in the transport equation, which is 
then simplified to: 

V • [iL(E;„f)V/e+] + ] + Q{E,,r} = 0 , (27) 

where we have approximated the kinetic energy of the positron by its total energy 
Ef.. Positrons with energies in the GeV-TeV range lose energy in their propagation 
through the galactic diffusive halo mainly via the inverse Gompton scattering (ICS) 
on the interstellar radiation field (ISRF) and via the synchrotron losses on the 
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Galactic S-field: b = feics + bsyn- The rate of energy loss due to ICS is given by 

bics{E,, E) = I de 2 dE^ {E^ - e) " V iSRF(e, E) . (28) 

Here da^^/dEj denotes the differential cross section of inverse Compton scattering 
of a positron with energy Eg, where an ISRF photon with energy e is up-scattered 
to energies between E^ and E^ dEj. In the limit e, me ^ Eg, kinematics and 
the neglect of down-scattering require that e < E^ < {1/Ee + l/dygc)”^ = 
with 7 e = Eejme- Besides, /isRF(e,fO is the number density of photons of the in¬ 
terstellar radiation field, which includes the cosmic microwave background, thermal 
dust radiation and starlight. An explicit model of the interstellar radiation field 
can be found, e.g., in Ref. For a positron energy of E^. = I GeV, 5ics ranges 
between 4.1 x 10“^^ GeVs“^ and 1.9 x 10“^^ GeVs“^, depending on the position in 
the Galaxy. On the other hand, the rate of energy loss due to synchrotron emission 
is given by 

6sy„(Ae,r) = ^crT7e^ , (29) 

where ctt = 0.67 barn denotes the Compton scattering cross section in the Thomson 
limit and B^/2 is the energy density of the Galactic magnetic field, which is con¬ 
ventionally taken as R ~ 6p,G exp(—12:|/2 kpc — r/10 kpc),l^ although the size and 
spatial dependence are not precisely known. At the position of the Sun this magnetic 
field yields a synchrotron loss rate given by 6syn — 4.0 x 10“^^(Re/ GeV)^ GeV s“^. 
In many analyses, the rate of energy loss is approximated by a spatially constant 
function parametrized by b[E) = with Eq = 1 GeV and te = 10 ^® s. 

Most experiments measure the positron fraction, defined as the flux of positrons 
divided by the total flux of electrons plus positrons, which is less susceptible to 
systematics since most sources of systematic error, such as detector acceptance or 
trigger efficiency, cancel out when computing the ratio of particle fluxes. This is the 
case of CAPRICE,!^ HEAT,I^ AMS-01,1^ PAMEL^lllll and AMS-0211I1I111 Some 
experiments also measure the positron flux itself, such as HEAT,I^ PAMELA^^ 
and AMS-02.1^ The most striking feature in the data is the rise in the positron 
fraction at energies larger than ~ 10 GeV, first observed by PAMELA an 
recently confirmed by the AMS-02 collaborationJHM^ While dark matter annihi¬ 
lations or decays provide a possible explanation for this observation, the cosmic 
positron backgrounds are still poorly understood and it is not possible at the mo¬ 
ment to make a definite statement about the origin of this excess. Nevertheless, 
the exquisite quality of the positron data can be used to set stringent limits on 
the dark matter annihilation cross sectionHUIlll Fig.[^ shows the limits on the an¬ 
nihilation cross section for the annihilation XX e^e~x which follow from the 
non-observation of a sharp feature in the positron fraction, and which approxi¬ 
mately read crn < 3 x 10“^®(3 x 10“^®)cm^s“^ for = 10(100) GeV.^^ The 
plot also shows for comparison the limits on the same annihilation channel from 
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DM coupling to electrons 



[GeV] 


Fig. 8. 95% C.L. upper bo und s on the annihilation cross section for the higher order processes 

XX e"^e “7 derived in Ref. |98| using the AMS-02 data on the positron fraction.^^ In the figures, 
it was assumed the Einasto profile and a dark matter particle quasi-degenerate in mass with the 
scalar mediator. For comparison, we also show as dashed lines the limits on the cross sections for 
the corresponding channels from gamma-ray data from the Fermi- LATJ^ and H.E.S.S. data.I^ 


the non-observation of a sharp gamma ray feature in the Fermi-LAT® and H.E.S.S. 
data.l^ As apparent from the plot, the current limits from gamma-rays are stronger 
than those that follow from the positron fraction data. 


4.4. Neutrinos 

Dark matter particles traversing the Sun could scatter-off a nucleus in the solar 
interior and lose energy. After several scatterings, the dark matter particles eventu¬ 
ally sink to the solar core where they accumulate.^i^SI The subsequent annihilation 
of the dark matter particles produces, depending on the final state, a high energy 
neutrino flux in the direction of the center of the Sun which could be detected in a 
neutrino telescopeli^ldSZJ Cosmic ray interactions with the solar corona constitute 
an irreducible background in this search, however, the predicted flux at the Earth is 
far below the sensitivity of current neutrino telescopes,therefore, the observation 
of an excess of high energy neutrino events correlated to the direction of the Sun 
would constitute a strong hint for dark matter annihilations. 

The time evolution of the number of dark matter particles N in the solar core 
is described by the following differential equationP^ 

^=Tc-CaN^-CeN, (30) 

at 


where Fc is the capture rate, Ca the annihilation constant and Ce the evaporation 
constant. For dark matter masses above ~ 10 GeV the evaporation of dark matter 
particles can be safely neglectedElllH^ Then, after solving Eq.(30), one finds an 
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annihilation rate as a function of time given by 

TA{t) = ^CAN{tf = ^Tc tanh^ (t/r) , (31) 

where 

The annihilation rate reaches a maximum when t ^ r. In this regime, captures 
and annihilations are in equilibrium and the annihilation rate reaches its maximum 
value, which is determined only by the capture rate: Ta = Tq/2. For thermally 
produced dark matter particles, this regime is never reached in the interior of the 
EarthP^ Instead, the annihilation rate is very suppressed due to the small number of 
dark matter particles trapped inside the Earth and, correspondingly, the expected 
high energy neutrino flux. Therefore, the detection of a signal from dark matter 
annihilations in the direction of the center of the Earth becomes very challenging, if 
not impossible in practice. On the other hand, for dark matter particles which have 
a sizable coupling to quarks, equilibration is generically reached inside the Sun and 
hence the high energy neutrino flux is not suppressed, thus making the center of 
the Sun a prime target to search for dark matter annihilations. 

In the class of scenarios under consideration the annihilation channels with 
largest branching fraction are yy —)■ ffV an yy —)• VV, with V, V' gauge bosons, 
which produce in their decay and hadronization a neutrino flux. It is important 
to note that the annihilations take place in dense medium of the solar core, there¬ 
fore the neutrino spectrum generated can significantly differ from the corresponding 
spectrum in vacuum, due to the interactions with the medium of the charged lep¬ 
tons and the hadrons before they decay. Concretely, the muons and light hadrons, 
such as pions and kaons, are stopped in the Sun before decaying and hence produce 
neutrinos with energies below 1 GeV. On the other hand, taus and heavy hadrons, 
such as charmed or beauty hadrons, decay in flight after losing a fraction of their 
kinetic energy, hence producing energetic neutrinos, while weak gauge bosons de¬ 
cay promptly producing neutrinos which typically dominate the energy spectrum at 
the highest energies. Finally, the neutrinos produced in the annihilation propagate 
from the solar interior to the surface, undergoing flavor oscillations and scatterings 
off solar matter, and from the solar surface to a neutrino telescope at the Earth, 
undergoing flavor oscillations in vacuum. 

The non-observation at IceCube of a significant neutrino excess in the direction 
of the Sun with respect to the expected atmospheric background then allows to set 
constraints on the capture cross section in scenarios where the dark matter particle 
couples to a light fermion EHEHIinil The resulting limits on the spin-dependent (left 
panel) and spin-independent (right panel) cross section are shown in Fig. for 
dark matter annihilations mediated by couplings to right-handed electrons (upper 
plots) and right-handed up- or down-quarks (lower plots) as a function of the dark 
matter mass, for eight different values of the parameter ranging between 
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Aimihilation into leptons: limits on crgQ 



Aiinihilatioii into light (juarks: limits on ergo 



Annihilation into leptons: limits on agi 



m,|OeVl 


Aiiniliilatioii into light quarks: limits on ergi 



Fig. 9. 90% C.L. limits on the spin-dependent (left plots) and spin-independent (right 

plots) interaction cross sections as a function of the dark matter mass for couplings to 
the right-handed electron (top plots) or to the right-handed first generation quark (lower 
plots), for various values of the mass ratio m^lniy^. The best limits on the scattering cross 
section from direct detection experiments are also shown for comparison. 


1.01 and 10 , where the dominant annihilation channel is /rJrZ and 7 Z respec¬ 
tively. The figure also shows the limits on the spin-dependent interaction cross 
sections from direct detection experiments COUPF^^^ and SIMPLE,*i 2 SI and on the 
spin-independent interaction cross sections from XENONlOCp^ and LUX.*ii4<'or 
couplings to leptons, the limits on the spin-dependent cross section are significantly 
stronger than the direct detection limits, while for quarks they are weaker, since 
in this case the dark matter particle annihilates mostly into final states involv¬ 
ing gluons which have a much larger branching fraction than those involving weak 
gauge bosons. On the other hand, the limits on the spin-independent cross section 
from IceCube are much weaker, at least one order of magnitude than the direct 
detection limits for couplings to cr and at least two orders of magnitude for cou¬ 
plings to Ur. For lighter dark matter particles, complementary constraints on the 
cross-section can be obtained from the non-observation at Super-Kamiokande of the 
neutrino flux produced by the stopped muons and charged pions and which peaks 
at energies below 1 GeVilililliHl 

A search for a high energy neutrino flux from dark matter annihilations has also 
been conducted in other regions of the Universe with an overdensity of dark matter 
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particles, such as the Milky Way center,lii2lhJll dwarf galaxiesthe Andromeda 
galaxy,IIi^ the Coma Clustei^^i^ and the diffuse extragalactic background lilHlliJo] 
For annihilations into W~^W~, the upper limit on the cross-section derived from 
neutrino data is 0(10“^^) cm^ s-i,*^^^ which is weaker than the corresponding limits 
derived from gamma-ray or antiproton data, unless the dark matter mass is in the 
multi-TeV range. Similar conclusions are expected for the annihilation channels into 
"fZ and ffZ, and which are the dominant ones in this scenario. 


5. Direct detection 



Fig. 10. Feynman diagrams for the processes contributing to the scattering of the dark matter 
particle off nucleons for a colored mediator (upper row) and an uncolored mediator (lower row). 
Diagrams with flipped initial and final state x are not shown for brevity. 


The rapid experimental development in direct detection makes it possible to 
probe even scenarios where the expected scattering rate of dark matter on nuclei is 
relatively small. It turns out that this is indeed the case for the simplified models of 
Majorana dark matter considered here. Nevertheless, one should keep in mind that 
the constraints depend on the local density and velocity distribution of dark matter 
particles, and therefore are subject to large uncertainties which are complementary 
to those affecting indirect constraints. In the following we assume a standard dark 
matter halo model characterized by a truncated Maxwellian distribution with vq = 
220km/s, Vesc = 544km/s and local dark matter abundance po = 0.4GeV/cm^ (see 


e.g. Refs. 20 109). 


5.1. Colored mediator 

The formalism for computing direct detection event rates is well developed and 


we do not repeat it here (see e.g. Refs. 121 122). Instead we briefly discuss the 
effective WIMP-nucleon interaction for Majorana dark matter which couples to the 
light Standard Model quarks via a scalar mediator with purely chiral coupling, and 
collect the quantities entering the cross section. The relevant Feynman diagrams 
are shown in Fig.|10[ For spin-independent interactions, the leading dimension six 
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operators Xl^XQlgQ both vanish, due to the chiral interaction and the 

Majorana condition, respectively. Consequently, the leading contribution arises from 
dimension eight operators,^^^ 

,2 

;{xi^D^x){qRlgD^qii - {D^qn)XgqR) ■ (33) 


/•SI 

■^eff-tree 


y 


2{m'i — (to^ + mqYY 


'^Tj V''^X ^ ''^QJ 

Additionally, a loop-induced interaction with the gluon content gives rise to an 
effective dimension seven interaction.^^ For spin-dependent scattering, the axial 
dimension six operator gives the dominant contribution, 

,,2 


/’SD 

^eff ,tree 


y 


^xTibxqRigibqR 
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Note that the tree-level contributions are enhancement for small mass splitting. The 
above expressions are applicable for niq — 3> rrip (see Appendix B). 


The corresponding scattering cross sections on nucleons are given by 
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where /xat is the WIMP-nucleon reduced mass, Aq^ is the spin-content of the nu¬ 
cleon ,^1^ and 


tto = 


y 


1 


(36) 


8 rriq — {m^ + ntqY ’ 

is the coefficient of the axial-vector effective interaction. The effective nucleon cou¬ 
pling for spin-independent scattering can be parametrized as 

= — Y X/ (/t^9 + 39 ( 2 ) + 3g(2)) g, - -I--G(2)gG , (37) 

q—u,d,s 


where cx {N\qq\N) and f^Q = 1 — d s fxq related to the integrated 

quark and gluon content,G(2), q(2), q{2) are the second moments of the gluon, 
quark and antiquark distribution,^!^ respectively, and 

%P' 1 

9<!~~~^7ZA2 TZf ’ (^^) 


8 (to2 _ (to^ _|_ m,)2)2 ’ 

is the coefficient of the effective dimension eight operator contributing to spin- 
independent scattering. The loop-induced interactions with the gluon content are 
proportional to the loop functions 

as 




as y , 2 T \ 
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(40) 
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where the loop integrals In are given iijiSSl (^ffQ ^se the corrected expressions from 


Ref. 124). As was emphasized in Ref. 124, the loop-induced effective couplings are 
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regular for rriri —> rn^ +mq, while the effective couplings and gq which arise from 
tree-level exchange of the mediator are resonantly enhanced. 

Due to the absence of dimension six operators for spin-independent scattering, 
limits on spin-dependent scattering yield the dominant constraint over a wide region 
of the parameter space. Nevertheless, since limits on spin-independent scattering are 
much stronger, the latter are important especially for small mass splitting, due to 
the less severe suppression of the dimension eight operators. 

In our analysis we consider spin-independent limits from XENONlOCp^ and 
from which are most sensitive at present within the mass range rriy, > 

40 GeV considered here. Although xenon is not an optimal material for spin- 
dependent searches, the presence of the isotopes ^^®Xe and ^^^Xe together with 
the large exposure and small background allow to set stringent limits on the spin- 
dependent WIMP-neutron scattering cross section, e.g. based on data collected by 
XENONlOOliHIIliUlFor the simplified models considered here, the coupling strength 
to neutrons and protons is comparable. Therefore, also experiments probing the 
spin-dependent proton cross section {e.g. SIMPLE,^^ COUPpE^ yield compara¬ 
ble constraints to those derived from the XENONIOO limits. Additionally, these in¬ 


teractions are probed by observations of the neutrino flux from the sun (cf. Sec. 4.41. 

As an example, we show in Fig. El the spin-independent as well as spin- 
dependent scattering cross section that corresponds to a thermally produced dark 
matter particle for a scenario where the mediator couples to right-handed up-quarks. 
The cross section depends strongly on the mass splitting and also on the dark matter 
mass. For small mass splitting rriq/m^ = 1.1(1.01), the cross section is suppressed 
below < 200 GeV (ITeV) due to strong coannihilations, while being sizable 
for larger masses due to the resonant enhancement of the scattering cross section. 
For 0(1) mass splitting, the qualitative behavior changes due to the absence of 
coannihilations. In addition, spin-dependent scattering is more sensitive compared 
to spin-independent scattering, since the latter is induced by dimension eight oper¬ 
ators for Majorana dark matter. For the standard halo model, thermally produced 
dark matter can be excluded from constraints on spin-independent scattering from 

LU}(!nni 

for 230GeV< < 550 GeV if rriq/m^ = 1.1, and from spin-dependent 

scattering constraints from XENONIOCP^ for < 150 GeV if = 2. 

The expected increase in sensitivity in XENONIT,*^^ and later LUX-ZEPLIN 
(LZ), will allow us to probe a significant fraction of the accessible parameter space 
for thermally produced Majorana dark matter that couples to first-generation 
quarks. If dark matter couples to heavy quarks, the constraints are significantly 
weaker, e.g. about an order of magnitude in ctsi cx; for bottom quarks.^^^ Nev¬ 
ertheless, even if dark matter couples only to top quarks, LZ and XENONIT are 
sensitive to scattering cross sections expected for thermal production for mass split¬ 
tings in the fewxC>(10%) range.^ 
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DM coupling to u-quark (spin independent) 



DM coupling to u-quark (spin dependent proton) 


DM coupling to u-quark (spin dependent neutron) 




Fig. 11. 90%C.L. direct detection limits for spin-independent scattering (upper plot), spin- 

dependent scattering off protons (lower left plot), and off neutrons (lower right plot) for 
mrj/mx = 1.01,1.1,2,5. Black lines correspond to the scattering cross section expe cted 
for t herm al production, and red l ines to experimental constraints from XENO N 100, 
LUX,E5I COUPF^SIl and SIMPLE.E^ Also shown is a projection for XENON ITliHHHlUl 


5 . 2 . Uncolored mediator 

It is interesting to note that even for a pure coupling to leptons, dark matter inter¬ 
actions can be potentially probed in next-generation direct detection experiments 
due to the loop-induced nucleon couplingUlliHlElll In particular, the charged medi¬ 
ator gives rise to an effective electromagnetic coupling. However, for Majorana dark 
matter, the loop-induced dipole moments vanish, and the leading contribution is 
the electromagnetic anapole moment. 


Leff = , 


(41) 
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DM coupling to leptons (anapole moment) 



[GeV] 


Fig. 12. 90%C.L. direct detection limits on the anapole moment A/ hm derived from LUX 

limits.i^ Black lines correspond to the anapole moment expected for thermal production 
for = 1.01,1.1,2,5, respectively, if dark matter couples to right-handed muons. 

Also shown is a projection for XENONllEHElIland for LUX-ZEPLIN.i^ 


witlP 


y^e f 3 / r^\ _ ml + 3w^ - 3mj 
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(42) 


This expression is valid if the momentum transfer is smaller than the mass of the 
particles in the loop, in particular for £ = r. For £ = e, a form factor depending 
on the momentum should be included.i^ 

Although the anapole moment is enhanced logarithmically for very small mass 
splitting, it is difficult to probe thermally produced dark matter in this region due to 
efficient coannihilations, which largely reduce the value of the coupling y required for 
thermal production. For a coupling to muons, the anapole moment corresponding 
to thermal production is shown in Fig. |12| t ogether with constraints inferred from 
LUXl^as well as prospects for XENON ITllHHlIIlIl and for LZ.^^ Taking coannihilations 
into account, the expected sensitivity of LZ may be large enough to detect thermally 
produced leptophilic dark matter in a small portion of the parameter space. The 
LUX constraints are significantly weaker compared to Fermi LAT and FI.E.S.S. 
limits except for very fine mass splittings ~ 1.01, where they are comparable 

(c/. Fig.[^. Prospects for XENONIT are competitive with GAMMA-400 and CTA 
for similarly small mass splitting and low dark matter masses. 
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6. Collider signatures 

The non-observation at LEP or the LHC of signals of new physics allows to set limits 
on the parameter of the dark matter model we consider in this review. The search 
strategies crucially depend on the mass difference between the dark matter particle 
and the scalar mediator. When the mass difference is very large, only the dark matter 
particle can be directly produced at the collider experiment. In this case, the most 
relevant experimental signature is a monojet, monophoton or mono-W/Z boson plus 
missing transverse momentum, which is generated by the pair production of dark 
matter particles which recoil against a jet or gauge boson. This signature has been 
studied in detail by ATLAdUlIHl and CMS,E^ mostly within an effective operator 
descriptiorEmEHl (see Refs |l40| - [143| for some recent discussions of the validity of 
this approach). On the other hand, when the mass difference is moderate, the scalar 
mediators can be directly pair-produced IllESlllHEllllililtll®! The subsequent decay of 
the mediator into the dark matter particle and the corresponding quark or lepton 
leads to signatures with missing transverse momentum and two jets, or with two 
same-flavor, opposite-charged leptons, respectively. Signatures of this type are being 
studied extensively in connection to SUSY searches |iS2H153] Lastly, when the mass 
difference is small, the scalar mediator is also directly produced in the collider. 
However, the jets or charged leptons produced in their decay are too soft to be 
detected llMHlSfil therefore the final state is completely invisible to the detector and 
the mono jet or monophoton constraints give again the strongest limits on the model 
parameters. 

The production rate of mediators at the LHC, and correspondingly the limits on 
the model, crucially depends on whether they carry color charge or not. We discuss 
in what follows each case separately. 

6.1. Colored mediator 

We consider first the case of a colored scalar mediator which couples to the dark 
matter particle and to the right-handed up or down quark. The Feynman diagrams 
contributing to the production at the LHC of the mediator are shown in Fig.[T3| 
The first five diagrams are mediated by strong interactions and contribute to the 
subprocesses gg —>■ r^ry, for the first four, and qq —>■ rjfj, for the fifth. The sixth and 
seventh diagrams are important when the dark matter particle couples to a first 
generation quark, and correspond to the subprocesses qq —)■ gfj and qq —>■ 7777 via the 
exchange of a dark matter particle in the t-channel (note that the latter process is a 
consequence of the Majorana nature of our dark matter candidate). Finally, the last 
two diagrams correspond to the subprocess qq —>■ 7777 via electroweak interactions 
and give a negligible contribution at the LHC for a colored mediator. The produc¬ 
tion cross sections for the different subprocesses can be found in Fig.[^ for some 
specific choices of the parameters and -y/s = 8 TeV. It follows from the plot that for 
small Yukawa couplings the dominant production process at the LHC is gg rjfj, 
while for larger Yukawas (77 > 0.5 for the parameters of the figure) the dominant 
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Fig. 13. Feynman diagrams for the processes contributing to the direct production of a colored 
mediator and an uncolored mediator (lower row) at the LHC. 


process is uu —>■ 7777 , due to the enhancement of the rate by the parton distribu¬ 
tion functions. This process is specific for a Majorana dark matter particle, and the 
corresponding cross section scales with the squared of the dark matter Majorana 


mass (c/. right plot in Fig. 141. Thermal dark matter production typically requires 


a sizable Yukawa coupling, cf. Fig. therefore for this particularly interesting case 
the dominant production mechanism of colored scalar mediators at the LHC is due 
to the exchange of a dark matter particle in the t-channel. This is in contrast to the 
squark production in the Minimal Supersymmetric Standard Model, which is dom¬ 
inated by channels involving the strong interaction. As a consequence, the ATLAS 
and CMS limits on simplified SUSY scenarios cannot be straightforwardly applied 
to our toy model, since the efficiency of the search depends on the underlying hard 
process. A critical discussion can be found in Ref. 

A dedicated search for the signals of Majorana dark matter coupling to a light 
quark has been performed in Ref. from a re-analysis of the ATLAS searcfP^ 
for jets and missing energy based on £ = 20.3fb“^ at -y/s = 8 TeV (see also 


Refs. 144 -147 for related works). The resulting upper limits on the pair-production 
cross section of the colored mediator are shown in Fig.[^ as a function of the mass 
splitting — 1 and various dark matter masses. For small mass splitting, the 

constraint reaches a plateau given by cr ~ 0.6 — 30 pb, depending on the dark mat¬ 
ter mass TOy = 800 — 200 GeV. For 0{1) mass splitting, the constraints are more 
stringent, of the order 1 — 0.01 pb. The theoretical expectation for a thermally pro¬ 
duced dark matter particle is also shown in Fig.fl^ The coupling strength required 
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Fig. 14. Production cross section of pairs of the colored mediator at the LHC as a function of 
the Yukawa coupling y (left plot) and the dark matter mass (right plot), for different production 
channels. 

for thermal production can be excluded for example for 1 < — 1 7 if 

Tn^ = 200 GeV and for 1 < mr^jm^ — 1 < 2 if = 500 GeV. 

These exclusion limits are obtained under the assumption that r; —is the 
only possible decay channel, which is the case as long as the simplified model is 
a good description. When relaxing this assumption, and considering an additional 
hypothetical untagged decay channel, the multijet limits degrade due to two effects: 
first, due to the reduced event rate suppressed by BR(7y —>• x<z)^- Second, for very 
large total width, the narrow width approximation becomes questionable.*^^ 

Qualitatively, this picture remains the same when considering multiple colored 
mediators and/or a coupling to down-type quarks, or left-handed first generation 
quarks. For example, for the case of two mediators that couple to the right-handed 
up- and charm-quarks, respectively, the quantitative change in the excluded Yukawa 
coupling strength is less than about 30% for m,-^ > SOOGeV or > 2, but can 

be larger otherwise due to the two times larger contribution from 5't7(3)c-mediated 
processes to the production cross section.*^ 

Apart from searches for multiple jets, constraints on a colored t-channel mediator 
have also been derived for the monojet channel.li^ShliZl addition to the process 
99 xX9i '^19' radiation of the gluon off the initial state or off the mediator in 
the t-channel, also the processes qg —> XX 9 , via initial-state radiation or on-shell 
production of xVj contribute to the monojet channej^ For very small mass splitting, 
the direct pair production of the mediator also contributes to the monojet signal.*^^ 
For Majorana dark matter with 0(1) mass splitting of the mediator, the constraints 
on the Yukawa coupling y inferred from monojet searches have been found to be less 


'^When radiating add itional gluons, these processes can also give a certain contribution to multijet 
searches (c/. Ref.lSHfor an analysis for Dirac dark matter). However, for Majorana dark matter, 
their relative contribution is less important due to the large contribution from qq rjrj. 
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= 200 GeV = 300 GeV 







= 800 GeV 



Fig. 15. 95%C.L. constraints on the production cross section for direct production of the 

colored media tor f or = 200,300, 500, 800 GeV, obtained from a reinterpretation of the 
ATLAS searclPil for jets and missing energy. Blue solid lines show the upper limit when 
using two hard jets in the matching procedure, and blue dashed lines when using only 
one jet. The black dotted line corresponds to the theoretical prediction for y — 0, and the 
black solid line for y = The shaded regions correspond to an estimate of uncertainties 
(see Ref. 1^ for details), and the blue dots are the upper limits obtained by ATLAS for a 
simplihed supersymmetric model. 


stringent compared to the multijet channel EStliZi pgj. degenerate masses, the 

monojet exclusion derived in Ref. |154[ taking rifj production via strong interactions 
into account, covers some parameter space which is complementary to the multijet 
search discussed above, which is however already excluded by XENONlOO.I^^Similar 
constraints can be derived also from monophoton searches .1^^ 


6.2. Uncolored mediator 

For a mediator which is singlet under SU{3)c, the production at LHC is clearly 
much less efficient than for colored mediators, and therefore the corresponding con¬ 
straints are expected to be significantly weaker. A scalar mediator which couples the 
Majorana dark matter particle to a Standard Model lepton can be pair-produced 
via the Drell-Yan process (c/. diagrams in the lowest row of Fig. 13), and the sub- 
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sequent decay into a pair of opposite-sign, same-flavor leptons as well as missing 
transverse momentum carried away by the dark matter particles provides a clear 
signature. Since the production processes as well as the decay is in this case identi¬ 
cal to the simplified supersymmetric model containing neutralino and sleptons, the 
corresponding exclusion limits apply directly in this case, and no re-interpretation 
is necessary. On the other hand, this also implies that collider constraints depend 
only on the masses and and are essentially independent of the coupling y 
(as long as y is large enough such that the mediator decays promptly, which is the 
case possibly except for m,, — < Wr, see e.g. Ref. 157). 

The constraints obtained from the LHC are sensitive to > lOOGeV, 

because smaller mass splittings would lead to leptons which are too soft to be effi¬ 
ciently discriminated from backgrounds. The production cross section depends on 
whether the mediator couples to right- or left-handed leptons, due to the different 
coupling strength to the Z boson. Searches for direct slepton production have been 
performed both by the ATLAS and CMS collaborations ll^HHi^For example, the AT¬ 
LAS search yields a rather mild constraint ranging up to > 250(300)GeV and 

^ 100(150)GeV in the right-(left-)handed case, for > lOOGeV. Nev¬ 
ertheless, these constraints are complementary to LEP limits, > lOOGeV, 

which apply for smaller mass splitting. The current constraints are summarized in 

Fig.[T9l 


7. Complementarity of searches 

Each strategy to search for traces of non-gravitational interactions of dark mat¬ 
ter, via cosmic rays from annihilation, scattering off nuclei, and production of dark 
matter particles at high-energy colliders, provides pieces of information that are 
complementary in two different respects: first, their sensitivity depends on different 
properties of the dark matter particle and its interactions with the Standard Model 
particles. Second, each channel is affected by different types of systematic uncer¬ 
tainties. These include above all the dark matter density and velocity distribution, 
both locally at the Earth as well as for the primary targets of indirect searches, in¬ 
cluding the Galactic Genter and dwarf galaxies as prominent examples. In addition, 
the ability to separate a signal from dark matter annihilation from astrophysical 
foregrounds and efficiently reject backgrounds in direct detection represent major 
challenges. The systematic uncertainties can be reduced, for example, by improved 
dynamical constraintspSE^ progress in the theory of structure formation, as well 
as by exploiting the different spatial and spectral morphology of astrophysical fore¬ 
grounds with respect to a dark matter annihilation signal. The spectral feature 
from internal bremsstrahlung is a particular example for a characteristic signature 
for which uncertainties from astrophysical foregrounds are considerably reduced. 

In the following, we discuss the complementarity of indirect, direct and collider 
searches for the simplified Majorana dark matter models subject of this review. 
To make a quantitative comparison, we assume the standard halo model for direct 
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[GeV] 


[GeV] 


Fig. 16. Compilation of current constraints on Majorana dark matter coupling to the right- 
handed up-quark, expressed in terms of the Yukawa coupling y, for various values of the mass 
splitting ATLAS constraints on jets and missing energy are shown in green, as well as 

direct detection (XENONIOO in red, LUX red dashed) and indirect detection (Fermi, H.E.S.S. in 
blue). The red dotted lines shows prospects for XENONIT, blue dotted for GAMMA-400 and blue 
dashed for CTA, respectively. The thick black lines corresponds to thermal production. 


detection as described in Sec.[^ and quantify gamma-ray constraints relative to the 
Einasto profile from Sec.|^ As discussed before, it is important to keep in mind the 
different sources of systematic uncertainties in each case. 


7.1. Complementarity of constraints on the coupling 

In Fig. |16[ we show a compilation of constraints and prospects discussed in the 
previous sections, for the case of a colored mediator. All cross sections have been 
converted into limits on the Yukawa coupling strength y, and the thick black line 
shows the coupling for which thermal freeze-out yields a relic abun¬ 
dance ~ 0.12. For very small mass splitting = 1.01 direct detection 

constraints from LUX and XENONIOO are by far dominant due to the resonant 
enhancement of the scattering cross section in this limit. LUX limits exclude ther¬ 
mally produced dark matter for a small range of masses between 1 — 2 TeV, and 
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m, |GeV] [GeV] 


Fig. 17. Same as Fig.|16| but for coupling to the right-handed muon. 


XENONIT will probe deeply into the region of thermal production. For mass split¬ 
ting niri/m^ = 1.1, direct detection constraints are much weaker, but still dom¬ 
inant for my_ ^ 1 TeV. For higher masses, gamma ray constraints from internal 
bremsstrahlung are formally the best constraint, although at unrealistically high 
values of the Yukawa coupling. Constraints from jets and missing energy at LHC 
are comparable to LUX constraints around « 0.5 — 1 TeV. For even larger mass 
splitting rrirf/m^ = 2, the overall picture changes and collider constraints become 
dominant over the whole range of masses > 100 GeV. This is due to two ef¬ 
fects: first, there is no resonant enhancement of the direct detection cross section 
in this regime. Second, the larger mass splitting leads to more energetic jets in the 
decay rj ^ XQ which make the discrimination against background more efficient. 
The limits inferred from the ATLAS search for jets and missing energy are slightly 
stronger than the coupling required for thermal production for « 200 — 600 GeV. 
Constraints from internal bremsstrahlung also become weaker, because this process 
is suppressed for large mass splitting. For very large splitting = 5, collider 

constraints are dominant as well, and exclude thermally produced dark matter for 
TOx ^ 300 GeV. 
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The constraints shown in Fig. 16 correspond to a coupling to right-handed up- 
type quarks. Nevertheless, the qualitative picture remains the same for a coupling 
to quarks of the first two generations HHH For example, for a coupling to down- 
type quarks the direct detection constraints remain essentially unaffected, the col¬ 
lider constraints weaken only slightly due to the dependence of the production 
cross section on the parton distribution in the proton, and limits from internal 
bremsstrahlung weaken by a factor due to the smaller electric charge of down- 
type quarks. Similarly, an extension to a model containing two mediators which 
couple to the up- and charm-quark, respectively, leaves collider constraints essen¬ 
tially unaffected, because the dominant production process for Majorana dark mat¬ 
ter is sensitive to the valence quarks.^ Note however that the presence of the second 
mediator influences thermal freeze-out and therefore affects the coupling for which 
the density matches the Planck value.^^It is also important to note that collider and 
direct detection constraints are considerably weaker for a coupling to the third gen¬ 
eration, while indirect constraints are identical in the limit ^ ^t/b 


For an uncolored mediator which couples dark matter to leptons, constraints 
from indirect detection are dominant for > 100 GeV at present. Therefore, 
searches for a spectral feature from internal bremsstrahlung represent an important 
probe of this class of models. The corresponding constraints and prospects are 
summarized in Fig.|17| 


7.2. Complementarity for thermal production 


When assuming that the Majorana particle x constitutes the dominant form of 
cold dark matter, and that it was produced via thermal freeze-out, one parameter 
of the model can be fixed by requiring = 0.12. Using this constraint to fix 

the Yukawa coupling strength, y = yth(rn^,mrf), the parameters are completely 
specified by the dark matter and the mediator mass, respectively. As discussed in 
Sec.[^ requiring thermal production as well as perturbative couplings implies that 
only a finite region in this parameter space is theoretically viable. In the following 
we discuss the current experimental constraints, and to which extent the viable 
region will be covered in the near future. 

For a colored mediator, collider searches as well as direct detection constraints 
exclude dark matter masses up to = 600GeV, see Fig. 18 (upper panel). How¬ 


ever, the exclusion regions are almost fully complementary: constraints from jets 
and missing energy are relevant for mass splitting > 2, and those from 

direct detection for smaller splittings. As discussed previously, this is due to an 
interplay of a resonant enhancement of the scattering cross section off nuclei on 
the one hand, and the efficiency of producing jets with large transverse momentum 
on the other hand. We also show monojet constraints which are sensitive to the 
quasi-degenerate region.ll^ Gamma ray constraints from internal bremsstrahlung 
at present are sensitive to couplings larger than those required for thermal produc¬ 
tion, when assuming an Einasto profile. The blue contour lines shown in Fig.|18| 
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show the ratio of the excluded annihilation cross section to the thermal cross sec¬ 
tion. The regions inside the contour lines are excluded if the gamma ray flux from 
annihilation was enhanced relative to the Einasto profile by a corresponding factor, 
either due to a cuspier prohle or due to substructures. 

The lower panel of Fig. 18 shows prospects for XENONIT and CTA. The former 
experiment will be able to probe a significant fraction of the theoretically viable 
parameter space in this model. Depending on the control of systematic uncertainties, 
CTA will be significantly more sensitive than H.E.S.S.. Furthermore, the region of 
highest sensitivity, in the multi-TeV region and for small mass splitting of order of 
tens of percent, will not be covered by direct or collider searches. Therefore, CTA 
will provide complementary information, although a detection may only be expected 
for optimistic assumptions on the dark matter distribution in this model. 

For the case of an uncolored mediator, constraints and prospects are shown in 
Fig. 19 Collider constraints from LEP and LHC partly exclude the region < 
100 GeV. The blue contour lines correspond to gamma ray constraints similarly 
as before. The reason why the sensitivity peaks at a mass splitting of 0(10%) is 
due to an interplay of two effects: first, the strength of the internal bremsstrahlung 
feature increases with smaller splitting. Second, the thermal coupling yth(rn^,mr^) 
decreases for m^, due to coannihilations. Although dark matter couples only 

to leptons at tree-level in this scenario, it may be possible to probe a small region of 
the parameter space with the future direct detection experiment LUX-ZEPLIN, due 
to the loop-induced anapole moment.^ Nevertheless, indirect detection will remain 
the most sensitive probe for large parts of the parameter space. 


8. Conclusions 

We have reviewed the phenomenology of a model where the dark matter particle is 
a Majorana fermion that couples to a light Standard Model fermion via a Yukawa 
coupling with a scalar mediator. The simplest scenario contains only three free pa¬ 
rameters, the dark matter mass, the mediator mass and the Yukawa coupling, and 
leads, for appropriate choices of the parameters, to the observed dark matter abun¬ 
dance via thermal freeze-out. The small number of free parameters of the model 
allows us to systematically scan the parameter space, explore its rich phenomenol¬ 
ogy, and investigate the complementarity among the various search strategies. The 
most characteristic aspects of the model, and which make its phenomenology rather 
unique, are: 

• Requiring dark matter thermal production and perturbativity of the cou¬ 
plings leads to a finite parameter space, that could potentially be closed by 
experiments. 

• The annihilation process yy —>■ // controls the dark matter abundance, 
however, it has a negligible cross section at present times. Instead, the 
dominant annihilation channels today are the higher order two-to-three 
process yy —>■ ffV and the one loop process yy —>■ VV', with V, V gauge 
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DM coupling to u-quark 



DM coupling to u-quark (prospects) 



Fig. 18. Compilation of current constraints from collider searches, direct detection (XENONIOO, 
LUX) and indirect detection (H.E.S.S.) on thermally produced Majorana dark matter with colored 
mediator, coupling to right-handed up quarks. The lower panel shows prospects for direct and 
indirect constraints achievable with XENONIT and CTA, respectively. 
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DM coupling to leptons 



[GeV] 


DM coupling to leptons (prospects) 



Fig. 19. Compilation of current constraints and prospects for thermally produced Majorana dark 
matter with uncolored mediator, coupling to right-handed muons. 
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bosons. 

• When the dark matter particle and the scalar mediator are very degenerate 
in mass, the annihilation process XX ffl leads to a sharp gamma-ray 
spectral feature which resembles a distorted gamma-ray line, and which 
could be detected in gamma-ray telescopes. The detection of such sharp 
feature in the gamma-ray sky would unambiguously point to a dark matter 
origin. 

• The tree level annihilation processes xx —>■ ffZ and the one loop processes 
XX -A- ZZ and xX 7-^ necessarily have a non-vanishing cross section 
if the dark matter particle couples to a Standard Model fermion. These 
processes lead to an exotic contribution to the antiproton flux, even for 
leptophilic dark matter particles. 

• For dark matter particles coupling to light quarks, the scattering rate with 
nucleons is enhanced in the mass degenerate limit, due to the resonant 
s-channel exchange of the scalar mediator. Away from the resonance, the 
spin-independent rate is highly suppressed since the coefficients of the di¬ 
mension six operators vanish. Therefore, in this regime, the spin dependent 
scattering ought not be neglected. In fact, current experiments probing the 
spin-dependent interaction yield the dominant constraint over a wide region 
of the parameter space. 

• For dark matter particles coupling to quarks, the colored mediator plays 
a pivotal role in colliders, since it could be light enough to be directly 
produced in parton collisions. For a Majorana dark matter particle which 
couples to a light quark, the total production cross section can be largely 
enhanced due to the contribution from quark-quark collisions. After being 
produced, the scalar mediator decays into a dark matter particle and a 
quark, which produces a jet. When the mass difference between the dark 
matter particle and the scalar mediator is large, the final state consists of 
two or more jets plus missing energy. On the other hand, when the mass 
difference is small, the produced jets are too soft to be triggered and then 
the final state is invisible. However the process could be detected from 
the emission of a gluon off any of the colored states involved, producing a 
monojet signal. In either case, a correct description of the phenomenology 
of the model at colliders requires the inclusion of the colored mediator. 

Within the scenario of a colored mediator, the theoretically viable parameter 
space for a thermally produced dark matter particle will be probed to a large extent 
by next-generation experiments, while the scenario of an uncolored mediator is 
much more challenging. Indirect, direct and collider constraints are covering largely 
complementary regions of parameter space, depending both on the dark matter mass 
and the mediator mass. Consequently, a null result at the next LHC run would not 
preclude the possibility to observe a signal at XENONIT. Similarly, a negative 
result by XENONIT would leave open the possibility to observe a signal from 
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internal bremsstrahlung at CTA. Nevertheless, in an optimistic case, a signal could 
be observed in more than one channel if the mass splitting between the mediator and 
the dark matter is of order one, and a combination of the characteristic signatures 
mentioned above would allow to pinpoint this dark matter model. 
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Appendix A. Annihilation cross sections 


We include in this appendix the expressions for the relevant annihilation and coan¬ 
nihilation cross sections in the toy model discussed in Section]^ in the limit nif —t 0 
and keeping the lowest order in the expansion in the relative dark matter velocity 
V (unless otherwise specified). In the following formulas, y is the Yukawa coupling 
between the Majorana dark matter particle y, the scalar ry, and the right-handed 
Standard Model fermion /a, which has electric charge qf and color Nc- The cross- 
sections for annihilations into /r/r were calculated in Ref. 160 for 77 and gg via 
a one-loop diagram in Refs. |161||162| for 'yZ in Ref. 


for ZZ in Ref. 

Ir/rI in Refs. |^|2^ for /r/rZ in Refs. |^|^|7^ and for fRfRg in Ref. 
cross sections for coannihilations were for example derived in Ref. mi 
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The 


Two-to-two annihilation into fermions 


Ncy^Ji- 4mj/s / ^2 

avff = ---^ m) + —- ^[Sml - 

32 TT f -I- J \ 12 (_ ^2 _(_ ^2 j 

— 5rnf{m^ — mj)'^ + 2m^{4:m* — + 0(v^) | (A.l) 


where s = = 4 . Note that we did not expand the square root in the velocity 

to avoid a singular behaviour that occurs for the particular case 1 — rnj/m'^ < v^. 
If this is not the case, the square root can be expanded in the velocity as well. 
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Two-to-two annihilations into gauge bosons via loops 
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with Li2(a:) = J2T=i the dilogarithm function, while A-yz and Azz are defined 
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Co being a Passarino-Veltman function. These expressions satisfy 
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Two-to-three annihilations 
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The spectra of gauge bosons are obtained by integrating the differential cross-section 
over the fermion energy, with integration limits given by j- 

— My) 12. The total cross-section can be obtained by integrating over the 
remaining energy with limits = My and Ey’^^ = + My/{Am^). The 

corresponding expressions for annihilation into left-handed fermions can be found 
in Ref. EH 

For the case of annihilations into a fermion-antifermion pair and a photon, the 
total cross section is given by: 
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Coannihilations 
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Appendix B. Effective operators for direct detection 

We want to obtain the effective Lagrangian for the scattering with a parton, xq —)■ 
Xq, and eventually for the scattering with a nucleon, by integrating out the scalar. 
Its equation of motion reads 

+ ml)g = -yqRX, = -yx^R ■ (B.l) 


By integrating out the scalar at tree level, one obtains the effective interaction 
Lagrangian 

The momentum of the scalar in the scattering is given by ~ 

+ mq)^. Therefore, we would like to expand the denominator around 
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By iterating this relation, we can rewrite the expansion as 
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where acts on qji and on Finally we also use the Fierz identity in the 
form 

(Xte)(te)x = ^X1^XQR'l^iQR - ^X7^75X9fl7M75te , (B.7) 

where we have used that qa = PnqR is chiral. Putting everything together, and 
including the terms coming from the hermitean conjugate terms, one can formally 
write the effective Lagrangian as a sum of spin-independent (SI) and spin-dependent 
(SD) contributions. 
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The lowest orders correspond to the effective Lagrangian of Drees and Nojiri.l^^ 
For the spin-dependent one, the lowest order is 

(B.IO) 


y 


^eff = -^^X7'-75Xte7M75to: 


which corresponds to Eq (1) of 164 for a chiral interaction with a = —b = y. 
The n = 0 order of the spin-independent interaction vanishes because X7^X = 0 
for Majorana spinors (however, e.g. xi^^vX in non-zero). Using this property, the 
lowest non-zero contribution can be written as 

7/2 
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2 (Am 2)2 
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(B.ll) 


This corresponds to Eq (10) of 164, again for a chiral interaction, and includes the 
contribution from the twist-2 operator. The expansion parameter is schematically 
given by 




2m^mq 


Am 2 


Amy^{Eq + mq) 

Am 2 


(B.12) 


where Eq is a typical energy scale for the parton inside the nucleon (more precisely, 
the higher orders contain higher moments of the parton distribution function in the 
nucleon). When the scalar and the dark matter particle are nearly degenerate, the 
expansion parameter is thus of the order of m^Eq/(m'^ — (m^ + mq)^) ~ Eq/ (m,, — 

m^ — mq). 
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